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We report results of susceptibility x and ^Li NMR measurements on LiVSi206. The tempera- 
ture dependence of the magnetic susceptibility x{T) exhibits a broad maximum, typical for low- 
dimensional magnetic systems. Quantitatively it is in agreement with the expectation for an 5 = 1 
spin chain, represented by the structural arrangement of V ions. The NMR results indicate antifer- 
romagnetic ordering below Tn — 24 K. The intra- and interchain coupling J and Jp for LiVSi2 06, 
and also for its sister compounds LiVGe2 06, NaVSi2 06 and NaVGe2 06, are obtained via a modified 
random phase approximation which takes into account results of quantum Monte Carlo calculations. 
While Jp is almost constant across the series, J varies by a factor of 5, decreasing with increasing 
lattice constant along the chain direction. The comparison between experimental and theoretical 
susceptibility data suggests the presence of an easy-axis magnetic anisotropy, which explains the 
formation of an energy gap in the magnetic excitation spectrum below Tn, indicated by the variation 
of the NMR spin-lattice relaxation rate at T ^ Tn. 
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I. INTRODUCTION 

In compounds of the series AVX2O6, where A=Li,Na 
and X=Si,Ge, the trivalent V ions occupy the centers of 
VOg octahedra forming chain-like structured elements. 
Because of this arrangement of the S = 1 magnetic mo- 
ments , these materials may in principle be considered as 
a physical realization of an ensemble of one-dimensional 
spin S = 1 chains with an antiferromagnetic intrachain 
couphng J. These, according to Haldane [ij, are ex- 
pected to adopt a non magnetic ground state, separated 
by an energy gap from magnetically excited states. 

Nevertheless the measured magnetic susceptibility x 
of the first member of this family, LiVGe206 0, was 
found to be clearly different from the expected x(T) 
of spin S = 1 chains. The discrepancy was first at- 
tributed to some sophisticated types of intrachain in- 
teractions 0, 13 , but subsequent ^Li NMR experiments 
01 revealed a three-dimensional antiferromagnetic order 
below Tn — 24 K, most likely due to a small but non 
negligible coupling Jp between moments of neighbouring 
chains. Magnetic susceptibility and ^^Na NMR measure- 
ments later confirmed a similar behaviour of NaVGe2 06 
with Tn = 18 K Q. 

It is conceivable that the compounds of the series 
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AVX2O6 represent a physical realization of solids with 
spin 5 = 1 chains but with interchain interactions that 
are large enough to quench the expected Haldane gap 
and to provoke a magnetically ordered ground state, as 
was suggested in Ref. on the basis of magnetic sus- 
ceptibility and specific heat measurements. 

In this article, we present the results of measurements 
of the magnetic susceptibility and of the ^Li NMR re- 
sponse of LiVSi206. Inspecting the temperature depen- 
dence of the magnetic susceptibility reveals no anomaly 
that would indicate the onset of magnetic order. How- 
ever, an abrupt broadening of the NMR spectra upon 
cooling reveals a transition from the high-temperature 
paramagnetic to the low-temperature antiferromagneti- 
cally ordered phase at Tn — 24 K. In addition, a promi- 
nent peak in the temperature dependence of the NMR 
spin-lattice relaxation rate T^^(T) at the same temper- 
ature supports this interpretation. 

Including previously obtained experimental data we 
extend our discussion to cover the series AVX2 06- Us- 
ing a modified random phase approximation (RPA), we 
determine the values of J and Jp from the correspond- 
ing values of Tn and T,„ax, where the latter temperature 
corresponds to the maximum of the magnetic suscepti- 
bility. While Jp is approximately constant, J varies by 
a factor of 5 across the AVX2O6 series. We then com- 
pare the measured magnetic susceptibilities to those cal- 
culated by quantum Monte Carlo (QMC) methods with 
corresponding values of J and Jp, and argue that the 
observed discrepancies may be explained by a single ion 
anisotropy. 

The paper is organized as follows. In Section ^1 we 
address the crystal structure and the synthesis of the 
LiVSi206 sample. Section lTlll is devoted to presenting the 
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results of magnetic susceptibility measurements, and in 
Section lTVI we display the results of the NMR measurents. 
In Section we determine J and Jp for the compounds 
in the AVX2O6 series, and compare the experimental and 
calculated magnetic susceptibilities. In the Appendix we 
present the technical aspects of the QMC calculations 
used in this paper in some detail. 



II. CRYSTAL STRUCTURE AND SAMPLE 

The crystal structure of LiVSi2 06 is schematically 
shown in Fig. ^ Typical for the AVX2O6 series, we note 
the chains of isolated VOg octahedra joined at the edges. 
These chains are linked and kept apart by Si04 tetrahe- 
dra. Considering the most likely oxidation states of Li"*" , 
Si''+ and 0^~, the V ions are expected to be trivalent. 
The magnetic moments are thus due to two electrons oc- 
cupying the 3d shell of V'^"'" ions, thus forming an 5 = 1 
configuration. The effective interaction between these 
moments is mediated by the oxygen 2p-electrons, and is 
direct only for neighbours placed on the same chain. In 
this geometry, the intrachain nearest neigbour coupling J 
is expected to be much larger than any possible coupling 
Jp between moments on different chains. 

LiVSi206 crystallizes in a monoclinic structure, space 
group C2/c 0. The LiVSi2 06 powder sample was pre- 
pared with the following procedure. First, a precursor 
with nominal composition LiSi20x was synthesized by 
heating a mixture of Li2C03 (Aldrich, 99.99%) and Si02 
(Aldrich, 99.995%) at 725°C for 50 hours in air and a re- 
annealing at 730°C for 20 hours in air. Then a stoichio- 
metric mixture of LiSi20x and V2O3 (Aldrich,99.99%) 
was pressed into pellets and annealed at 990° C for 40 
hours in an evacuated silica tube. Subsequent heat treat- 
ment for 60 hours in another evacuated silica tube was 
necessary to obtain the desired composition of LiVSi2 06. 
According to x-ray powder diffraction data, the sam- 
ple contained 95% of the desired phase with the pyrox- 
ene structure and 5 mass percent of non-magnetic Si02. 
The lattice parameters of LiVSi2 06 were determined by 
room temperature X-ray diffraction measurements, and 
resulted in a = 9.634 A, 6 = 8.586 A, c = 5.304 A, and 
P = 109.69°. 

In Table m we list the structural parameters of the unit 
cell of the compounds in the AVX2O6 series. In Table Hll 
we also list the length of the intrachain V-V bonds, as 
well as the angle between two successive V-V bonds in 
the same chain. These values were obtained by using a V- 
V bond corresponding to the vector v — Saa + 6bb + 6cc, 
where a, b and c are the unit cell vectors. From the 
geometry of the chains it follows that Sa = 0, since the 
chains lie in a-planes, 6c — 0.5 and Sb — ±|(56|, where 
the sign switches because of the zig-zag character of the 
chains. The parameter \Sb\ is different for the individual 
compounds, and is calculated from the exact positions of 
the V atoms. We note only small variations of the V-V 
bond length and the V-V-V angle across the series. 




FIG. 1: Representation of the crystal structure of LiVSi2 06. 
The V^^ ions are located in the center of the dark-grey VOe 
octahedra, while the Si^^ ions occupy the centers of the light- 
grey Si04 tetrahedra. The Li^ions are represented by the 
spheres. 





Space group 


a (A) 


6(A) 


c(A) 




LiVSi206 




C2/c 


9.634 


8.586 


5.304 


109.69° 


NaVSi206 


m 


C2/c 


9.634 


8.741 


5.296 


109.90° 


LiVGeaOe 


\2] 


P2i/c 


9.863 


8.763 


5.409 


108.21° 


NaVGe2 06 


[5] 


P2i/c 


9.960 


8.844 


5.486 


106.50° 



TABLE I: Values of the unit cell parameters in the AVX2O6 
family. 



III. MAGNETIC SUSCEPTIBILITY 

Using a commercial SQUID magnetometer, we mea- 
sured the magnetization M of 18.8 mg of powdered 
LiVSi206 (number of mols N = 8.95 x 10^^), for tem- 
peratures T between 4 and 340 K and in magnetic fields 
fioH between 0.01 T and 5 T. Apart from a temperature 
independent component, the molar magnetic susceptibil- 
ity X — M/{HN) was found to be identical in different 
magnetic fields above 10 K. 

In Fig. |2]we display an example of x(2^)7 measured in 
a magnetic field fJ-oH = 5 T. The temperatures are too 
low to identify a high-temperature Curie- Weiss behaviour 
and thus to establish the oxidation state of the V ions. 
Reducing the temperature from 340 K, x(T) increases 
and reaches a broad maximum at T^ax = 115 ± 1 K. This 
type of temperature dependence of x is typical for low- 
dimensional magnetic systems. Upon cooling to below 15 
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5b (A) 


V-V length (A) 


V-V-V angle 


LiVSi2 06 


0.188 


3.105 


62.7° 


LiVGeaOe 


0.186 


3.158 


62.1° 


NaVGeaOe 


0.193 


3.256 


63.2° 



TABLE II: V-V bond data for selected compounds in the 
AVXaOe family. 
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FIG. 2: Magnetic susceptibility x of LiVSi206 as a function 
of temperature T, measured in a magnetic field H = 5T. The 
inset magnifies the temperature range 15-35 K. 



K, x{T) exhibits a smooth upturn, -which is attributed to 
the lo"w-temperature Curie- Weiss tail caused by a small 
amount of paramagnetic impurities. 

We note that in this data set there is no evidence for 
any anomalies in x(T) at temperatures belo-w Tmax- In 
particular, as is sho"wn in the inset of Fig. [21 x{T) varies 
smoothly bet"ween 20 K and 30 K, the temperature range 
in which, as reported in Sec. IIVI a magnetic transition 
is revealed by the NMR data. 

This is some-what surprising because for LiVGe206 |3, 
NaVGe2 06 and NaVSi2 06 a prominent kink in 
x(r) indicates the onset of a three-dimensionally antifer- 
romagnetically ordered state upon cooling. 



IV. 
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Li NMR 



The NMR experiments probed the same sample as that 
used for the magnetization measurements. The ^Li NMR 
spectra -were obtained by monitoring the integrated spin- 
echo intensity as a function of the irradiation frequency / 
in a fixed external magnetic field H. The selected values 
^iqH = 7.0495 T (7.05 T in the following) and ^iqH = 
4.5533 T (4.55 T in the following) were calibrated by 
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FIG. 3: ^Li NMR spectra of LiVSi206 for selected tempera- 
tures above 25 K. The integrated echo intensity, multiplied by 
the temperature T, is plotted as a function of the frequency 
shift K. In the upper graph, two spectra taken at T = 290 K 
but in different magnetic fields hqH = 4.55 T and hqH — 7.05 
T are compared. 



monitoring the resonance frequency in liquid D2O. 
The respective ^Li reference frequencies, /o = 116.643 
MHz and /o = 75.340 MHz, were obtained from /o = 
7/ioi?, with 7 = 16.5463 s~-'^T~^ as the gyromagnetic 
factor of the ^Li nuclei. The spin echo was generated 
with a two pulse 7r/2-delay-7r spin-echo sequence. 

Fig. 13 shows examples of recorded spectra for hqH = 
7.05 T and fioH = 4.55 T at selected temperatures above 
25 K. They were obtained by irradiating a frequency win- 
dow of about 20 kHz. The integrated echo intensity, mul- 
tiplied by the temperature T, is displayed as a function 
of the frequency shift 



K 



fo 



1 



(1) 



where / is the irradiation frequency. In the temperature 
range between 25 K and 294 K the shape and the total 
intensity of the NMR signal are essentially the same . 

For T — 293 K, we show two spectra measured in the 
two different magnetic fields of 7.05 T and 4.55 T. As 
shown in Fig|3 they coincide in their central parts but ex- 
hibit some differences, indicated by the arrows, at the up- 
per and lower ends of the signal. This is attributed to the 
first order quadrupolar splitting of the nuclear transitions 
between the Zeeman states m = ±3/2 and m — ±1/2 
0. This frequency splitting is field independent and, in 
the spectra in the higher field, is masked by the magnetic 
broadening of the NMR line, which is proportional to the 
applied magnetic field. The width in frequency space at 
the basis of the spectrum for ^qH = 4.55 T allows for an 
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FIG. 4: Relative frequency shift Kp^ak of the peak frequency 
/peak of the ^Li NMR line, displayed as a function of the 
magnetic susceptibility x- The data are for temperatures T > 
25 K and in magnetic fields H = 7.05 T and = 4.55 T. The 
dashed line represents the best linear fit to the experimental 
data. Inset: The dashed line again represents the best linear 
fit. 



estimate of the upper limit of the quadrupolar frequency 
/q < 80 kHz of the ^Li nuclei. 

Above 25 K the lineshape of the signal is compatible 
"with that of a randomly oriented powder with anisotropic 
shifts. The anisotropic part may be attributed to the 
dipolar coupling between the Li nuclei and the paramag- 
netic V-moments. Such lineshapes were already observed 
in the case of NaVGezOg and LiVGezOe Q. 

In Fig. ^ the y-axis corresponds to the shift ifpcak of 
the frequency /peak at the maximum of the NMR signal. 
For T > 25 K, -fCpcak is plotted as a function of x{T)- 
The data for the two different fields collapse onto the 
same curve, and -f^pcak(x) can reasonably well be fitted 

by 



-^peak — ^^pcakX ~^ '^pcak 



(2) 



indicating that the temperature dependence of the shift 
of the NMR signal is of purely magnetic origin. The 
hyperfine coupling corresponding to the signal peak is 



^poak = apeak • N^iB = 1280 ± 50 G 



(3) 



If the hyperfine coupling expected for a purely dipolar 
interaction of the V moments with the Li nuclei is calcu- 
lated as in Refs. |1,[3, the result is Af^^^ = 210 G. The 
discrepancy between this calculated and the observed 
value is attributed to an additional non-dipolar hyperfine 





^pcak (G) 


^0 (G) 


LiVSiaOe 
LiVGezOe [4] 
NaVGeaOs f5] 


1280 ± 50 

580 
1300 ± 50 


1050 ± 50 

480 
1140 ±50 



TABLE III: Values of ylpeak and Aq for selected compounds 
in the AVXaOs family. 



couphng Ao = 1070±50 G, such that Apeak = A"^'"^ 



peak 



-Ao. 



In Table ITTll we compare the values of Apeak and Aq for 
three compounds of the AVX2O6 series. 
In the inset of Fig. 2|we plot 



FWHM 



(4) 



where FWHM is the full width at half maximum of the 
NMR signal and /q the above mentioned reference fre- 
quency, again as a function of x- Above 30 K, the data 
can be accomodated by 



(5) 



The corresponding anisotropic part of the hyperfine cou- 
pling 



NfiB = 700 ± 100 G 



(6) 



is close to the value a'^^ — 



650 G if calculated in the 
same way as in Ref . 5j , assuming purely dipolar coupling 
between the V magnetic moments and the Li nuclei. 

The inset of Fig. ^also shows that the value of a at 25 
K is distinctly larger than predicted by Eq. (0), thus indi- 
cating the onset of antiferromagnetic correlations. For a 
discussion of the low-temperature behaviour we turn our 
attention to the NMR spectra at temperatures below 25 
K. Fig. displays examples of spectra measured at tem- 
peratures below 30 K and in a magnetic field (J-qH = 7.05 
T, recorded by irradiating a frequency window of about 
100 kHz. Above Tn = 24 K, a single NMR line charac- 
terizing the paramagnetic state is observed. Upon cool- 
ing its intensity decreases drastically between 25 K and 
24 K, and the narrow line has completely vanished at 
23 K. Below Tn = 24 K, a broad, rectangular shaped 
NMR signal develops with decreasing temperature. This 
shape is typical for ordered moments in a powder sample. 
The interpretation of an onset of magnetic order is also 
substantiated by the H- and T-dependence of the width 
A/ of the signal, which is represented as a function of 
temperature T in the inset of Fig. [S] for magnetic fields 
^iqH = 4.55 T and hqH = 7.05 T. Indeed, A/ is field 
independent and, for T ^ 0, we argue that A/ tends 
to saturate at a constant value of A/, = 3.4 ± 0.4 MHz. 
Thus, below Tn, the magnetic moments residing on the V 
ions adopt a three-dimensional antiferromagnetic order. 
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FIG. 5: ''Li NMR spectra of LiVSi206 for selected temper- 
atures below 25 K. The echo intensity is plotted as a func- 
tion of the frequency / for selected temperatures below 25 K. 
The measurement were performed in a fixed magnetic field 
H — 7.05 T. Inset: line width A/ as a function of T, in 
magnetic fields H = 4.55 T and = 7.05 T. 



The internal magnetic field Hiat at the Li sites attains a 
saturation value Hint = A/*/27 = 1030 ± 60 G. 

To gain further information on the phase transition at 
Tn = 24 K, we measured the temperature dependence 
of the spin-lattice relaxation rate, T^^{T). T^^ was de- 
termined by monitoring the recovery of the ^Li nuclear 
magnetization m after the application of a long comb of 
radiofrequency pulses. The experiments were performed 
in a magnetic field fioH = 7.05 T, irradiating a frequency 
window of about 100 kHz in the center of the NMR sig- 
nal. An exponential recovery 



m(t) =moo[l-exp(-VTi)] 



(7) 



with recovery time t was observed across the entire tem- 
perature range between 5 and 295 K. The single exponen- 
tial law is appropriate, since we simultaneously irradiate 
all the three possible transitions of the / — 3/2 Li nu- 
clear spins, with the maximum splitting of the order of 
/q « 80 kHz. In Fig. Elwe display the temperature de- 
pendence of the spin-lattice relaxation rate T^^. Above 
Tn, T^^ slowly decreases by a factor of three between 295 
K and 30 K upon cooling. A prominent peak in T^^{T) 
at Tn = 24 K reflects the magnetic phase transition. 

The temperature dependence of the spin-lattice relax- 
ation rate at low temperature provides additional infor- 
mation about the low-energy magnetic excitations. Be- 
low Tn, T]~^ decreases by orders of magnitude (a factor 
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FIG. 6; ■^Li spin-lattice relaxation rate Tj ^ of LiVSi2 06 as 
a function of temperature T. The maximum of Tj~^ at Tn 
is indicated by the dotted vertical line. Inset: ^Li spin-spin 
relaxation rate Tj"^ as a function of T. Tn is indicated by the 
dotted vertical line. 
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FIG. 7: ■^Li spin-lattice relaxation rate Tj^ ^ of LiVSi2 06 as a 
function of inverse temperature . The solid line represents 
the function T"^ = Cexp(-A/r) with A = 14.5 K. 



of 100 between 25 K and 10 K, as shown in Fig. [7|. indi- 
cating the formation of a gap in the spectrum of magnon 
excitations in the antiferromagnetically ordered phase. 
This interpretation is supported by Tj~^(T) well below 
Tn, which varies according to 



Tfi(T)^exp(-A/T) 



(8) 



as is emphasized by the solid line in Fig. [7| Because of 
the presence of a small amount of paramagnetic impuri- 
ties, the actual value of the gap may be somewhat larger 
than the estimated value A = 14.5 K. A comparison of 
gap values, all established from T^^{T) well below Tn, 
is made in Table HVl for three AVX2O6 compounds; they 
all turn out to be almost equal. Such gaps are usually 
attributed to an easy-axis single-ion magnetic anisotropy 
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A(K) 


LiVSiaOe 
LiVGeaOe [4] 
NaVGeaOe [5] 


14.5 
14.7 
12.5 



TABLE IV: Values of the gap in the magnetic excitation 
spectrum A, calculated from the low-temperature dependence 
of Tj~^, for selected compounds in the AVX2O6 family. 







T^-T^oo (10' 


LiVSi2 06 
LiVGe2 06 \£ 
NaVGeaOe [5] 


1.9 ±0.1 
2.2 ±0.2 
1.0 ±0.1 


1.4 ±0.2 
1.2 ±0.2 
0.6 ±0.1 



TABLE V: Values of T^^t-^o ^nd TaTr^oo ^r selected 
compounds in the AVX2O6 family. For the calculation 
of T2 in LiVGe206, Eq. Q was replaced with m{t) — 



^ol|, an aspect discussed in more detail in Sec. 

Finally, we consider the temperature dependence of 
the spin-spin relaxation rate Tj"^. The corresponding 
experiments were made under the same conditions as the 
Ti experiments and T2 was determined from fits to the 
decay of the echo intensity m with the exponential law 

m(i) = moexp(-<7r2) . (9) 

Here t — 2t and r is the delay between the tt/2- and 
TT-pulses which generate the spin echo. 

In the paramagnetic phase above Tn, ^2""'^ is approxi- 
mately constant at a value « (1.4 ± 0.2) x 10^ s~^. As 
is shown in the inset of Fig. the phase transition is 
revealed by a sudden reduction of T2^(T) at Tn upon 
cooling. This decrease is followed by an upturn at even 
lower temperatures. As T ^ 0, T^^ tends to saturate to- 
wards the constant value T, 



(1.9±0.1)xl0^ This 

particular temperature dependence of Tj"^ seems to be 
typical for the the AVX2O6 series (see the T-^iT) data 
in In Table IVI we display the values of T^^^q and 



and notice that the values are of the same order 
of magnitude. 

In magnetically ordered states, T2~^ is often caused 
by magnon mediated interactions between nuclei. The 
unexpected upturn in T2^{T) upon cooling below Tn is 
incompatible with this scenario because the large gaps 
listed in Table HVl would lead to a rapid decrease of T2^. 
No such decrease has been observed, in particular for 
LiVGe206 down to 1 K, as reported in Ref. Pj. Therefore 



it must be that other mechanisms are essential for the 
spin-spin relaxation. 



V. WEAKLY INTERACTING SPIN S = 1 
CHAINS IN (Li,Na)V(Ge,Si)206. 

The results of the present and former studies H, ^ |^ 
clearly indicate that in the (Li,Na)2V(Ge,Si)206 series, 
the onset of the low temperature Haldane phase is in- 
tercepted by the developement of three-dimensional an- 
tiferromagnetic order. It is known that even a very 
small interchain coupling may quench the Haldane gap 
and induce a three-dimensional antiferromagnetic order 
[Til IT^ . In order to establish the ratios Jp/J be- 
tween interchain and intrachain couling for the different 
compounds, we offer an analysis of the x(T') data invok- 
ing quantum Monte Carlo calculations based on a model 
Hamiltonian. 

We model the system of magnetic moments at the V 
sites by 5 = 1 spins placed on a cubic lattice. The cou- 
pling to the next neighbours along the chain direction 
is antiferromagnetic, J > 0. Each spin is allowed to in- 
teract with the four nearest spins of neighbouring chains 
(see Fig. assuming an interchain coupling constant 
Jp. The spin S — 1 operators of the i-th spin on the 

The Hamiltonian is the sum 



r-chain are denoted by Sr 
of two terms. 



(10) 



where the first term represents the intrachain Hamilto- 
nian for a collection of antiferromagnetic S = 1 chains, 



H 



V,z+1 (J > 0) 



(11) 



The second term in Eq. (|10|) is the interchain Hamilto- 
nian 



Hj^ = J, 



(r,s) 



(12) 



where (r, s) indicates a pair of neighbouring chains. 
In the random phase approximation (RPA) 



Jp/J 



1 



(13) 



where yf is the staggered susceptibility of an isolated 
S ~ 1 chain (Jp = 0), which was calculated with QMC 
methods and reported in Ref. (see also Appendix 
IXj) . The factor 4 in the denominator corresponds to the 
number of neigbourhing chains. The validity of Eq. I|13|) 
was shown in Ref. 13] for Jp/ J < 0.2 (see also Appendix 
provided that the renormalized value of ^ = 0.695, 
instead of the classically expected ^ = 1, is used. From 
x(T) as obtained by QMC calculations, Tmax/J may be 



7 



0.3 



0.2 



0.1 



0.0 



S=1, modified RPA 
S=1/2 

LiVGe206 



NaVGe206 



NaVSi206 



LiVSizOe 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Tn / ''"max 

FIG. 8: Jp/J as a function of Tj^/Tma.^, resulting from the 
modified RPA approximation. The value of TN/Tmax for the 
AVX2O6 compounds are indicated. The dotted line represents 
Jp/J(TN/rmax)for weakly interacting spin S = 1/2 chains as 



given in Ref. 

^^ax — 0. 



k) for 

m 



Note the upper limit of Jp/J ~ 0.02 for 





Tmax (K) 


Tn (K) 


Jp/J 


J/fcfl (K) 


Jp/fcs (K) 


LiVSiaOB 


115 ± 1 


24.0 ± 0.5 


0.034 


85.7 


2.91 


NaVSi206 


48 


17.4 


0.068 


35.2 


2.40 


LiVGeaOs 


62 ± 1 


25 ±0.5 


0.079 


45.3 


3.58 


NaVGeaOe 


25.0 ±0.5 


18.0 ±0.5 


0.20 


17.6 


3.54 



TABLE VI: List of magnetic parameters of the AVX2O6 fam- 
ily that are relevant in the comparison between calculations 
and experiments. 



expressed as a function of Jp/J, using the approximation 



Tlnax/J = a2(Jp/J) ±ai(Jp/J)+ao 



(14) 



with 02 = —0.95, ai = 0.69 and aq — 1.32 (see Appendix 
IXj) . Combining Eqs. H13|l and H14|) . the relation between 



Jp/J and TN/T,„ax 

may be established. It is represented in Fig. (Sjand allows 
to derive the values of Jp/J for the compounds of the 
series AVX2O6. The results arc listed in Table lVIl These 
ratios differ from those given in Ref. j6j because the latter 
were obtained using the formula given in Ref. 0| , which 
is valid for spin S = 1/2 chains, and is represented by the 
dotted hne in Fig. El 

While the value of the interchain coupling Jp is almost 
constant across the series, the value of the intrachain 
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FIG. 9: Magnetic susceptibility x as a function of tempera- 
ture T of LiVSiaOg, compared with the x{T) resulting from 
quantum Monte Carlo calculations. Inset: dx/ dT{T). Note 
that the noise in the experimental dx/dT{T) data, computed 
from x{T) represented in Fig. |5| does not allow to confirm or 
exclude the presence of a small anomaly like the one reported 
in Ref. il. 



coupling J decreases by a factor of 5 from LiVSi2 06 to 
NaVSi206, LiVGe206 and NaVGezOg. The reason for 
such a variation cannot be established from the data at 
our disposal. It does not appear to be explained by geo- 
metrical features of the V-chains alone, since the relevant 
parameters differ only by small amounts (0.5% for the V- 

V bond length and 1° for the V-V-V angle; sec Table lll)l . 
It may therefore be that the position of atoms other than 

V and the resulting electronic structure are relevant for 
the indirect mediation of the magnetic interaction be- 
tween neighbouring V ions. Inspecting the positions of 
the O ions around the V ions, it turns ou that J decreases 
with increasing distortion of the VOe-ochtahedra. 

Next, we compare the experimental xC^) data of 
LiVGe206 3] and NaVGe206 with those calculated 
with QMC methods, employing the values of Jp and J 
from Table IVTl This is done in Figs. El El and llll where 
X is plotted as a function of temperature T. Bullets show 
the experimental data, and open circles the calculated 
X^^'^{T), adapted to the experimental data by varying 
the parameters xo and pcs in 

.QMC(j.) ^ ^ (iVAvMlK'ff Jxit = kBT/J) 



X 



(15) 



In Eq. (|15|) x a-nd t are the dimensionless susceptibil- 
ity and temperature, respectively, discussed in Appendix 
IXI Their relations to physical quantities is exempli- 
fied in Table IVIII The factor (A'av/^b) accounts for 
the correct units of a magnetic susceptibility. The ad- 
ditive constant xo represents the offset in the experi- 
mental data due to diamagnetic contributions and p^jj 
is a multiplicative factor, causing the high-temperature 
magnetic susceptibihty (iVAvMB)Pcff'5'('5'+ l)/3fcBr to be 
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FIG. 10: Magnetic susceptibility x a^s a function of tempera- 
ture T of LiVGe206, compared with the x{T) resuhing from 
quantum Monte Carlo calculations. Inset: dx/ dT{T). 
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FIG. 11: Magnetic susceptibility x as a function of tempera- 
ture T of NaVGe206, compared with the x{T) resulting from 
quantum Monte Carlo calculations. Inset: dx/dT{T). 



of Curie- Weiss type. We note that the value of PcB of 
NaVGe206 established in this way is somewhat larger 
than the one which resulted from the reported fitting 
of a Curie- Weiss law to the experimental xC^) data at 
high temperatures The insets of Figs. IHI El a-nd 
im represent the derivatives with respect to tempera- 
ture, {dx/dT){T), and {dx^^^ /dT){T). The discrep- 
ancy between the two can only partially be attributed 
to the presence of paramagnetic impurities, which might 
reduce dx/dT at low temperatures. The prominent peak 
in dx/dT for LiVGeaOe and NaVGeaOe (less so for 
LiVSi2 06) is not accounted for by the simulations, in- 
dicating that the Hamiltonian in Eq. I|1(J|) is not able to 



fully account for the magnetic properties of the AVX2O6- 
series. 

A more realistic model would include a term 

Hd = -DY,SI,SI, , (16) 

in the Hamiltonian of Eq. 1)10(1 , which describes a single 
ion magnetic anisotropy. D > favorizes the alignement 
of the moments along the z-axis, while D < favorizes 
configurations with magnetic moments in the plane per- 
pendicular to the z-axis. As discussed in the Appendix, 
a term with D > is appropriate for the present case. 
Indeed, it is the natural choice in order to account for the 
formation of the gap A in the magnetic excitation spec- 
trum, and to reproduce the sharp peak in dx/dT . This 
possibility was already considered in Ref. Further 
speculations would require a more detailed analysis of 
weakly coupled anisotropic spin S — 1 chains, which we 
leave for future investigations. Here, we only remark that 
Eqs. ifT^ and ifTH) have to be adapted if Z? ^ 0, which 
would yield values of Jp and J different from those given 
in Table ED 



VI. SUMMARY, CONCLUSIONS AND 
OUTLOOK 

Our experimental investigations confirm that in the 
AVX2O6 compounds the interchain coupling Jp is always 
large enough to quench the Haldane gap and provoke the 
onset of a three-dimensional antifcrromagnctic order at 
Tn > 0. While the calculated Jp is almost constant across 
the AVX2O6 series, the intrachain coupling J varies by a 
factor of 5. This variation does not seem to be correlated 
to geometrical parameters of the V chains. We therefore 
suspect that the electronic orbitals that are involved in 
the mediation of the coupling between the V magnetic 
moments along the chains, are affected by the positions 
of the other atoms in the structure. In particular, the 
growing distortion of the VOe octahedra seems to reduce 
the value of J. Band structure calculations and the sub- 
sequent determination of the Wannier functions of the V 
and O atoms of the four different compounds have to be 
made in order to understand the details of the intrachain 
V-V coupling. 

The comparison of the experimental and theoretical 
temperature dependences of the magnetic susceptibility 
indicates that the magnetic moments at the V atoms are 
affected by magnetic easy axis single-ion anisotropy, con- 
sistent with the observation of a gap in the magnetic 
excitation spectra below Tn- The origin and the mag- 
nitude of this effect could also be clarified starting from 
the results of band structure calculations. 
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Dimensionless 


Physical 


quantity 


quantity 


1 


J 


t 


T = J/kB ■ t 


jp 


Jp = J ■ jp 


X 


X = J-(iVAvMl)x 


5 


A = J-<5 


d 


D = J -d 


h 


H = J -h 


c 


C = C- NAvkB 



0.6 



TABLE VII: Conversion from dimensionless to physical quan- 
tities. 
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FIG. 12: Interchain coupling jp as a function of tN, deter- 
mined using the modified RPA according to Eq. ljA5^ (see 
also p. Il3|'). for different values of d > 0. The red symbols 
represent points {tN,jp), where tN is the temperature of the 
maximum specific heat c for a system with interchain coupling 
jp. Inset: Cp(t) for jp — 0.1 and different values of d > 0, cal- 
culated with QMC methods; the critical temperature tN is 
indicated. 



APPENDIX A: QMC CALCULATIONS, 
MAGNETIC SUSCEPTIBILITY AND CRITICAL 
TEMPERATURE 

In this Appendix we focus on the QMC calculations, 
aiming at clarifying some points of the discussion of the 
magnetic susceptibility in the AVX2O6 scries in Sec. Ivl 
For the QMC computations, the intrachain coupling is set 
to unity, i.e., J — 1, and all quantities are dimensionless. 
We employ this convention throughout this Appendix. 
Physical quantities with correct dimensions can be gained 
by using the corresponding energy values of J, according 
to Table ETH 

The QMC calculations are based on the Hamiltonian 

h = hj+jj^+hd , (Al) 

where the three terms are the intrachain Hamiltonian 
(see also Eq. (fTT|l ') 

the interchain Hamiltonian (see also Eq. H12|l 'l 

hj^ — jp ^ ^ 5*^,1; • Ss^i , (A3) 

{r.s):i 

and the magnetic anisotropy hamiltonian (see also (|16|l 'l 

hd = -dY,S'r,^S^,, . (A4) 



For nonintcracting chains {jp ~ 0), we considered a one- 
dimensional array of up to 100 S — 1 spins. For weakly 
interacting chains we used a cubic 8 x 8 x 20-lattice of 
S = I spins, with the largest extension corresponding to 
the chain direction. The spin S = 1 operators of the i-th 
spin on chain r are denoted by Sr,i- With (r, s) we denote 
a pair of neighbouring chains. We employed a stochastic 
series expansion [isl fl6| QMC code based on the ALPS 
library U^^. 

In the RPA approximation the effect of the interchain 
interaction hj is treated as an external perturbation 
of a single chain by the n neighbouring chains. Stan- 
dard mean field arguments show that the temperature 
of the transition from the high-temperature param- 
agnetic to the low-temperature antiferromagnetically or- 
dered state is related to the interchain coupling according 
to 

'"^WK^^) • ^^^^ 

The staggered susceptibility xtz P^r spin of an isolated 
S = 1 chain, described by the hamiltonian hj + hd, was 
calculated using the QMC technique for various values 
< c? < 0.5. By analyzing the results of QMC calcula- 
tions for the Hamiltonian IjAip on a cubic lattice (n = 4) 
with c? = 0, it was shown in Ref. that, for jp < 0.2, 
the interchain coupling renormalizes the value of ^ in Eq. 
(jA5|) from unity down to ^ = 0.695. In Fig. ^Iwe dis- 
play the result for jp{tj^) also for the case d > ,2Q||. We 
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FIG. 13: Magnetic susceptibility x calculated using QMC, 
as a function of the temperature t, for different values of jp. 
For jp = the location of tmax is indicated. Inset: tmax as 
a function of jp. The solid line represents the best fitting 
second order polynomial. 

















- 




X d = -0.)0 ■ 


// / 


— d = 0.00 ■ 


d * 

' A y 

o' ^ 


d = 0.07 - 


=' / 


^ d-0.)0 




d = 0.50 



■D 



0.8 
0.6 
0.4 
0.2 
0.0 



f for d = 

0.1: 0.5, 

0.0 : 1 



f 



FIG. 14: Xzz(^) (upper plot) and dxzz/dt{t) (lower plot), 
determined with QMC calculations, for fixed jp = 0.1 and 
different values of d. The transition temperature tN, corre- 
sponding to the peak in c{t), is indicated. 



find good agreement to discrete pairs of values (jp, iw) re- 
sulting of QMC computations using the full Hamiltonian 
(jAlj) . Such calculations yield, among other quantities, 
the temperature dependence of the specific heat c. For 
a given value of jp, the temperature of the peak in c is 
identified as (see inset of Fig. . 

We first focus on the case d — 0. In Fig. ^]we show 
examples of the temperature dependence of the magnetic 
susceptibility x(T), calculated using QMC, for selected 
values of jp. The temperature tmax, corresponding to the 
maximum of is indicated. The inset of Fig. I13l displavs 
imax(jp) which we find to be well approximated by 



a2jp 



with 02 = —0.95, ai = 0.69 and aq 
In this way, one obtains 



1.32. 



Up) 



tmax ( Jp ) 



(A6) 



(A7) 



where the numerator of the right hand side is the inverse 
function of jp(tN) given in Eq. ljA5|l . and the denomina- 
tor is given in Eq. ljA6p . The inverse of the function in 
Eq. 



Jp(iN/imax) 



(A8) 



is plotted in Figure |H1 in the form ( Jp/J)(TN/T'max)- 

Finally we consider the case of d 7^ which introduces 
an anisotropy of the magnetic susceptibility such that 



Xzz 7^ Xxx = Xyy Averaging over all directions, the 
effective value for a powder sample is 



_ 1_ 2_ 

Xpowdcr — ~^\zz H" ~^\xx 



(A9) 



The temperature dependence of Xzz(0 has been com- 
puted by QMC and is represented in the upper panel 
of Fig. ^] for jp = 0.1 and selected values of d. Note 
that Eq. HA6|I . and thus also HA8|I . do not hold if d > 0, 
since for fixed jp, the maximum of Xzz(i) is shifted to- 
wards higher temperatures. Nevertheless, two important 
conclusions can be drawn from the behaviour of Xzz (t) : 
First, d > induces an anomaly in Xzzit) at ^n, as is 
shown in the lower panel of Fig. ^1 This kink will be 
refiected in Xpowdcr (see Eq. HA9|l '). Secondly, for d > 0, 
the rapid decay of Xzz to zero upon cooling below indi- 
cates the formation of an energy gap d above the three- 
dimensional antiferromagnetically ordered ground state 
(in agreement with Ref. |lOj). The low-temperature be- 
haviour is consistent with an exponential suppression. 



Xzz{t) ~ exp{~S/t) 



(AlO) 



indicative of a gap S in the system. We leave a detailed 
analyis of the dependence of J on d and jp to future in- 
vestigations. Finally, we note that for d < 0, the presence 
of a gap S is excluded and in this case, Xzz tends to a 
non-zero value for f — > 0. 
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